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Abstract 

We establish a fundamental relationship between the averaged density of states and the extinction 
mean free path of wave propagating in random media. From the principle of causality and the 
Kramers-Kronig relations, we show that both quantities are connected by dispersion relations and 
are constrained by a frequency sum rule. The results are valid under very general conditions and 
should be helpful in the analysis of measurements of wave transport through complex systems and 
in the design of randomly or periodically structured materials with specific transport properties. 
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Fundamental questions in coherent transport of electromagnetic, electronic or acoustic 

riQ n 

waves [1|, |2j as well as applications to imaging in complex media [3|] have made wave prop- 
agation in random media a central issue in physics. Randomly or periodically structured 
materials allow to design media or devices with unconventional properties. This includes 
photonic crystals [4,], metamaterials for electromagnetic or acoustic waves ^] , strongly 
correlated disordered systems j^, 9], or materials generating non-diffusive transport 10]. 

The density of states (DOS) and the extinction mean free path (MFP) are fundamental 
concepts in coherent wave transport. The DOS shapes many macroscopic transport proper- 
ties The local density of states (LDOS) drives the spontaneous emission of light jl2| . 
and is at the root of novel imaging techniques base on field correlations [l^ • Fluctuations in 
the DOS or LDOS characterize the transport regime speckle patterns [l^ or the local 



structure of a complex medium 



16|. The extinction MFP 4xt, defined by Cxi = 4"^ + 



1 _ I ^-1 

"abs 



with 4 and Cbs the scattering and absorption mean free paths, describes the attenuation of 
the averaged (or coherent) field. When absorption is negligible, the extinction MFP equals 
the scattering MFP. The latter is an important quantity since kig is a measure of the strength 
of scattering, k being the wavenumber in the medium. The spectral behavior of 4 and the 
DOS were put forward in early studies of Anderson localization of light [17]. Moreover, the 
transition to localization in three dimensions is expected when kig ^ 1, according to the 
loffe-Regel criterion 1|. 

In this Letter, we establish a fundamental relationship between the averaged DOS in a 
random medium and the extinction MFP. As a consequence of causality and the Kramers- 
Kronig relations, we show that both quantities are connected by dispersion relations and 
are constrained by a frequency sum rule. We focus the derivation on light propagation 
in scattering media, but the results are valid under very general conditions and should be 
applicable to any kind of waves. 

Consider a scattering medium made of scatterers randomly distributed in free space (or in 
an otherwise homogeneous background medium). The (dyadic) Green function (or electric- 
field susceptibility) describes the response at point r, and at a given frequency u, to a point 
electric-dipole source p located at point r' through the relation E(r) = fiQUj"^ G{r,r')p. In 
free-space, the Green function reads : 



GnfR) = PV 



47iR 3kl 



2 



where I is the unit tensor, ko = uj/c, with c the speed of hght in vacuum (or in the ho- 
mogeneous background medium), R = r — r', R 



R| and PV stands for principal value. 
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In a random medium. 



Its Fourier transform is given by Go(k) = [{k"^ — /eq)! — kk] 
after averaging over the positions of the scatterers and assuming statistical translational 
invariance, the averaged Green function (G(r — r')) obeys the Dyson equation 



(G(k)) = Go(k) + Go(k) I](k) (G(k)) 



(2) 



in which S(k) is the self energy (or mass operator) containing the sum of all multiply 



connected scattering events 



18|. From Eq. ([2]) and the expression of Go(k), the averaged 



Green function can be written : 

(G(k)) 



(P - kl)I - kk - S(k) 
A(k) kk 



Fl-fc2e±(k) A;2Pell^(k) 



(3) 



where A(k) = I — kk/fc^ is the transverse projection operator. In the last equality, we have 
identified the effective dielectric function : 

S(k) 



eeff(k) =1 + 



k'^ 



(4) 



and its transverse and longitudinal projections e^(k) and e|!g(k). This quantity drives the 
propagation of the averaged (or coherent) field in the random medium, and is in general a 
non-local and anisotropic response function. In practice, determining the effective dielectric 
function is a difficult problem, that can only be solved under some (sometimes severe) 
approximations s], [3]. In the present study, we do not need to refer to a specific model. 
Our arguments rely only on the existence of the effective medium for the propagation of the 
averaged field. The latter is a direct consequence of the Dyson equation. 

Under the following hypotheses, (z) the random medium is isotropic on average, and 
{a) only fields variations on scales larger than the size of the scatterers and the correlation 
distance between scatterers are accounted for, the dielectric function becomes a local and 



isotropic (scalar) function 



Q, Q, i.e. 



-eff 



(k) = e|lfj(k) = eeflfl. We assume that these 



conditions are satisfied in the following. The averaged Green function in direct space now 
reads : 

exp(iA;eff-R) 



(G(R)) = PV 



I 



'^eff 
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where kpa = ko, with Ucs = ^/e^ the effective (complex) refractive index of the random 



medium 



20|. 



From the principle of causality, one can derive the Kramers-Kronig (K-K) relations that 
connects the real and imaginary parts of the susceptibility of any linear material. Regard- 
ing the optical response, the K-K relations are usually written in terms of the dielectric 
function [21^. It can be shown that the refractive index in passive materials 22 1 is also a 
quantity that satisfies the K-K relations (this amounts to showing that the refractive index 
is an analytic function in the upper half-plane of the complex frequency plane) 23(]. In a 
homogeneous medium, the imaginary part of the refractive index is associated with absorp- 
tion (or reflection in the particular case of a non-absorbing electron gas below the plasma 
frequency). In a scattering medium and below the homogenization threshold, even in the 
absence of absorption, the effective index has a non-vanishing imaginary part that corre- 
sponds to extinction of the averaged field by scattering. The lost energy is redistributed in 
the field fluctuations, whose averaged square modulus is the diffuse intensity. In terms of 
the effective index, the K-K relations read : 



Im n^fiiuj) 



1 + -PV 



uj' Im ncs{uj') 



uj'- 



2uj 
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(6) 
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As we will see, the K-K relations lead to interesting relationships between the spectra of 
the averaged DOS and the extinction MFP. Moreover, we will show that sum rules can be 
deduced under very general conditions. 

In an inhomogeneous medium, the LDOS at point r is given by p(u;, r) = 
2co'/(7rc^) ImTr G(r, r), where Tr denotes the trace of a tensor ll|, l2J]. The LDOS counts 
the number of radiation states in the frequency range [uj^uo + dcu], weighted by their con- 
tribution at point r. Note that in the near field, in the presence of evanescent states (due, 
e.g., to electromagnetic surface modes), the expression of the LDOS needs to be modi- 



fied 2a]. This subtlety is particularly important for the calculation of thermal radiation 
energy densities 26|. In a macroscopic random scattering medium with statistical transla- 
tional invariance, i.e. far from any surface boundary, the averaged LDOS does not depend 
on position and equals the averaged DOS : 

2a; 



p{uj) 



TTC" 



ImTr(G(R = 0)). 



(8) 



This quantity is the density of states that we consider in the present work. In the particular 
case of a vacuum, the DOS po{uj) is deduced from Go using the identity hm^^^o Im Go(R) = 
A;o/(67r)I One recovers the well-known result po{uj) = uo"^ / {ii'^ c^) . From Eqs. ^ and 

([H]), we obtain the averaged DOS in the random medium following the same approach, 
using the identity lim/j^o Ini (G(R)) = A;efr/(67r) I, which yields piu) = po{u!)Rencs{uj). 
The averaged DOS is given by the real part of the effective refractive index. We point out 
that this result has been obtained under the conditions of existence of an isotropic and 
local effective refractive index. This means that the effects of microscopic length scales on 
the order of the correlation length or the size of the scatterers have been disregarded. In 
particular, this expression coincides with the LDOS used to describe macroscopically the 
spontaneous decay of dipole emitters in homogeneous dielectrics. It is known that in dense 
materials local field corrections have to be incorporated in order to account for interactions 
on microscopic length scales 28|]. In the following, we shall refer to the DOS given above as 
the averaged macroscopic DOS. 

The second important quantity in our discussion is the extinction MFP £ext, that describes 
the attenuation of the averaged field by scattering and absorption [l| . More precisely, iext is 
defined as the decay length of the intensity of the averaged field. Therefore the extinction 
MFP is given by the imaginary part of the effective refractive index jiol . I20I : iext{^) = 
c/ [2ujlm n^fiiuj)]. 

The principle of causality implies a close connection between the averaged macroscopic 
DOS and the extinction MFP. Indeed, by a direct application of the K-K relations Eqs. ([6]) 
and (I7j), and using the expressions of p{uj) and iexti^^) given above, we obtain : 

f!^ i^lpvTll^,^' (9) 

Po{uj) vr Jq lu'^ - LO^ 

1 --i.^c^pv r "':?:"°''t? d^-. (10) 



These dispersion relations are the first result of this Letter. They are valid under very general 
conditions : passive and causal medium, and assumption of an isotropic and local effective 
medium for the description of the averaged field. They demonstrate that the averaged 
macroscopic DOS and the extinction MFP are not independent. Equation ([9]) shows that 
from the spectrum of the MFP, one can deduce the averaged macroscopic DOS (and vice 
versa using Eq. ffTOj) ). In practice, this means that from an extinction spectrum (a natural 
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measurement in spectroscopy), one could deduce the spectrum of the averaged macroscopic 
DOS. 

From the K-K relations, sum rules for the dielectric constant and the refractive index can 
be obtained {29;]. In particular, it is well established that the refractive index of any passive 
and causal medium satisfies ulmnesi^) [^GUesi^) — 1] dc<j = and f^lRenes^uj) — 
1] dcu = 29(1 . Beyond the principle of causality, the derivation of these sum rules relies on 



the assumption of a material behaving as a free-electron system in the high frequency limit : 
ee{i{uj) ~ 1 — ijJ^jijJ^ when — 00, where uj^ is an effective plasma frequency. Note that 
this hypothesis is not too restrictive, this high frequency behavior being expected as soon 
as the frequency is much larger than the resonant frequencies of the effective medium. The 
sum rules for the effective refractive index can be translated into new sum rules involving 
the averaged macroscopic DOS and the extinction MFP. In particular, we obtain : 

This relation is the main result of this Letter. It demonstrates that the spectra of the aver- 
aged macroscopic DOS and of the extinction MFP are intimately connected, and constrained 
by a simple sum rule. The simplicity and the generality of this relation are striking. Let 
us remind the three conditions of validity : (1) the medium has to be passive and causal, 
(2) the effective medium is described by an isotropic and causal dielectric function, (3) the 
medium behaves as a free electron gas in the high frequency limit. The second sum rule 
for the refractive index leads to a relation involving the averaged macroscopic DOS only : 
Jq°°[p(u;) — p<^{iS)\I uP' dcj = 0. This second sum rule was established previously in the context 
of spontaneous emission in dielectric media [s^. Regarding wave propagation in random 
media, it establishes a constraint on the potential modifications of the averaged macroscopic 
DOS. In particular, it shows that p(ll') is necessarily lower than the DOS in free space 
Po{uj) in a spectral range, and greater than po{uj) in another spectral range (the numerator 
p{uj) — po{uj) has to change sign somewhere for the integral to vanish). This behavior is also 
dictated by Eq. (fTT!) because ^xtl"^) > in a passive medium. In the following, we illustrate 
the general behavior induced by relations (I MTTI) in a particular case. 

Let us consider a random scattering medium with an effective dielectric function exhibit- 
ing a resonance at a particular frequency ujq. The resonance can be induced by an internal 
resonance of the scatterers, or of purely geometric origin (or both). We choose a Lorentz 
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model of the form : 

eefr(^) = 1 + , ^f. ^ (12) 
where the parameter T is an effective oscillator strength and F is the linewidth. From this 
expression, the effective index = is readily obtained numerically, as well as the 



spectra of the averaged DOS p{y>) and the extinction MFP £ext(t^)- As noted in Ref. 
a regime of strong scattering is identified when Imnefr(c<j) > Renefr(i-i-'), which corresponds 
to Hext < 1/2, where k = lie k^s is the wavenumber in the medium. In this regime, the 
effective medium satisfies Reeefr(ci;) < and the averaged field is strongly damped (the 
effective medium has a metallic character). In the regime kiext > 1/2, one has Iieees{uj) > 
and the effective medium has a dielectric character. The transition between these two 
regimes is driven by the parameter V = J-'ujq/T. For V < 2, one has lie ees{uj) > at all 
frequencies. For V > 2, there is a frequency range for which ReecfT(co') < 0, or equivalently 

k(^ext < 1/2. 

We show in Fig. [1] the spectra of the averaged macroscopic DOS (left) and of the extinc- 
tion MFP (right), in the case of an effective medium of dielectric character {V = 1). We 
observe that due to the resonance, the extinction MFP exhibits a minimum. Due to the 
K-K relationship, the averaged macroscopic DOS shows a typical dispersion behavior. This 
illustrate the fact that the DOS has necessarily a complicated structure around a minimum 
of the extinction MFP. In Fig. [21 we plot the same quantities for a medium with a strong 
metallic character {V = 20). The extinction MFP exhibits a region with a low value, corre- 
sponding to the spectral region for which ReeefT(to') < or kicxt < 1/2. In the same spectral 
range, the DOS is minimum, showing the existence of a pseudo-gap. Moreover, close to the 
lower pseudo-gap band edge, the DOS exhibits a strong oscillation that is a feature of the 
underlying dispersion relation. 

Finally, we note that in the case where the effective dielectric constant is known analyt- 
ically, one has Imegfr = 2 Im Re rzefr at each frequency. In terms of averaged DOS and 
MFP, this leads to the explicit relation lmees{^) = c/[ujiext{^^)] pi^^)/ Po{^)- This relation 
is useful when an analytical model of eefr is available, a situation that occurs only in a few 
specific cases. 

In summary, from the principle of causality, we have established dispersion relations and 
a frequency sum rule that connect the averaged macroscopic DOS and the extinction MFP 
in a random scattering medium. These relations constrain the spectral variations of both 
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FIG. 1: Normalized DOS p{uj) / po{io) (left) and extinction MFP l^xti^) in /^m (right) versus the 
normalized detuning from resonance 5 /T = (uj — loq) /T . The effective dielectric function eeff(w) is 
given by Eq. 1^, with loq = 10^^ Hz, F = 10^ Hz and V = I. 




6/r 6/T 
FIG. 2: Same as Fig. 1 with V = 20. 



quantities, and sustain general features of their spectral behavior close to a resonance of the 
effective medium. These results should be helpful in the analysis of measurements of wave 
transport through complex systems and in the design of randomly or periodically structured 
materials with specific transport properties. 
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